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Density matrix renormalization group (DMRG) is applied to a (l+l)-dimensional Ai^* model to study sponta- 
neous breakdown of discrete Z2 symmetry numerically. We obtain the critical coupling {\/ = 59.89 ± 0.01 
and the critical exponent j3 — 0.1264 ± 0.0073, which are consistent with the Monte Carlo and the exact results, 
respectively. The results are based on extrapolation to the continuum limit with lattice sizes L = 250, 500, and 
1000. We show that the lattice size L = 500 is sufficiently close to the the limit L ^ 00 pp. 



1. Introduction 

Hamiltonian diagonalization is a useful method 
for nonperturbative analysis of many-body quan- 
tum systems. If Hamiltonian is diagonalized, the 
system can be analyzed nonperturbatively at the 
amplitude level ^2JSA^ . However, in general quan- 
tum field theories, the method does not work 
without reducing degrees of freedom because the 
dimension of Hamiltonian increases exponentially 
as the system size becomes large. We need to 
find a way to create a small number of optimum 
basis states [5IB| . S. White proposed a power- 
full method called density matrix renormaliza- 
tion group (DMRG) PH). In DMRG, calcula- 
tion accuracy of target states can be controlled 
systematically using density matrices. White cal- 
culated very accurately energy and wavefunctions 
of Heisenberg chains composed of more than 100 
sites using a standard workstation. The cal- 
culation reproduced the exact value of ground- 
state energy in five digits or higher. DMRG has 
been applied to various one- dimensional models, 
such as Kondo, Hubbard, and t-J chain mod- 
els, and achieved great success. In many cases, 
DMRG can give more accurate results than quan- 
tum Monte Carlo. A two-dimensional Hubbard 
model has also been studied with DMRG in both 
real- and momentum-space representation |9I1U| . 
DMRG works well on small two-dimensional lat- 
tices and new techniques have been proposed for 



solving larger lattices DMRG has also been 
extended to finite-temperature chain models us- 
ing the transfer-matrix technique ■ In particle 
physics, the massive Schwinger model has been 
studied using DMRG to confirm the well-known 
Coleman's picture of 'half- asymptotic' particles 
at a background field 9 = n ^^l- It would be 
interesting to seek a possibility of applying the 
method to QCD in order to study color confine- 
ment and spontaneous chiral symmetry breaking 
based on QCD vacuum wavefunctions. 

DMRG was originally proposed as a method 
for spin and fermion chain models. In fermionic 
lattice models, the number of particles contained 
in each site is limited because of the Pauli prin- 
ciple. On the other hand, in bosonic lattice mod- 
els, each site can contain infinite number of par- 
ticles in principle. It is not evident whether 
Hilbert space can be described appropriately with 
a finite set of basis states in bosonic models. 
This point becomes crucial when DMRG is ap- 
plied to gauge theories because gauge particles 
are bosons. Before working in lattice gauge the- 
ories like Kogut-Susskind Hamiltonian ^1], we 
need to test DMRG in a simple bosonic model 
and recognize how many basis states are neces- 
sary for each site to reproduce accurate results. 
In this work, we apply DMRG to a At/)^ model 
with (1-1-1) space-time dimensions. We construct 
a Hamiltonian model on a spatial lattice. The 
model has spontaneous breakdown of discrete Z2 



1 



2 



symmetry and the exact values of the critical ex- 
ponents are known. We are going to justify the 
relevance of DMRG truncation of Hilbert space 
in the bosonic model by comparing our numeri- 
cal results with the Monte Carlo and the exact 
results |15I16| . Our calculations are done with 
lattice sizes L = 250, 500 and 1000. L = 1000 is 
about twice of the latest Monte Carlo one (15) . 




2. DMRG in the A0f+i model 

In the X(f)f^-^ model, we divide a Hamiltonian 

L L-l 

H ^^hn+^hn,n+l, (1) 

n—1 n—1 

into two parts 

The field operator 7r„ = a0„ is conjugate to 
[0m (0 1 '^nit)] = i^mn- The derivative has been re- 
placed with a naive difference. (Errors associated 
with the difference can be discussed if necessary 
\n\-) We rewrite Hamiltonian Q using creation 
and annihilation operators ajj and a„. 

1 i 
(i)„ = {al + a„) , 7r„ = — (aj^ - a„) , (2) 

where [am,ajj] = Smn and a„|0) = 0. Note that 
ajj and a„ are not creation and annihilation oper- 
ators in Fock representation. The index n of the 
operators aj^ and a„ stands for the discretized 
spatial coordinate, not momentum. Real-space 
representation is better for our purpose because 
local interactions are useful for DMRG. 

The finite system algorithm of DMRG is ap- 
plied to the Hamiltonian. A superblock Hamilto- 
nian Hs is composed of two blocks and one site: 

^^S = Hi^ + hn-l,n + hn + + Hr, (3) 

where Hj^ and H^i are effective Hamiltonian for 
the left and right blocks, respectively. 
(hn^n+i) is an interaction between the left (right) 
block and the inserted n-th bare site. The target 
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Figure 1. (A//i^)c is plotted as a function of 1/L 
for L = 250,500, and 1000. Extrapolation to the 
limit L^oo gives iX/^i'^)c = 59.89 ± 0.01. 



state is expanded as 

M N M 

i*)-EEE*»^-fcN'^'^)- (4) 

i=i j=i k=i 

where the index j is for the inserted bare site 
and i and k are for the renormalized blocks. 
The relevance of truncation with M and A'' 
is checked numerically by seeing convergence 
of energy and wavefunction. Parameter val- 
ues {M,N) = (10,10), which give good conver- 
gence, are used in all calculations. In Fig. 
and |51 we extrapolate the results to the con- 
tinuum limit and obtain the critical coupling 
{X/fi^)c = 59.89 ± 0.01 and the critical exponent 
(3 = 0.1264 ± 0.0073, which are consistent with 
the Monte Carlo iX/fJ.'^)c = 61.56t°;g and the 
exact (3 — 0.125 results, respectively. 

3. Conclusion 

We have determined the critical coupling con- 
stant (A//Lt^)c and the critical exponent f3 of the 
model by extrapolating the numerical results for 
finite but sufficiently large lattices to the contin- 
uum limit. DMRG truncation works well also in 
the bosonic model. The lattice with L = 500 can 
give results sufficiently close to the limit L —^ oo. 

The numerical calculations were carried on 
RIKEN VPP and Yukawa Institute SX5. 
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Table 1 



Various results for the critical 


coupling constant (A//i^ 


c are listed. 


Method 


Result 


Reference 


DMRG 


59.89 ±0.01 


This work 


Monte Carlo 


61.56lg1« 




Gaussian effective potential 


61.266 


in 


Gaussian effective potential 


61.632 


EH 


Connected Green function 


58.704 




Coupled cluster expansion 


22.8 < (A//i2)c < 51.6 




Non-Gaussian variational 


41.28 


ED 


Discretized light cone 


43.896, 33.000 


m 


Discretized light cone 


42.948, 46.26 
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Figure 2. is plotted as a function of 1/L for 
L = 250, 500, and 1000. Extrapolation to the 
limit L ^ oo gives /3 = 0.1264 ± 0.0073. 
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